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PROCEDE DE DETEI^IIia&TION DE L 1 IMPEDANCE D'UH SYSTEMS 

ELECTROCHH-HQUE 

Pour la determination des proprietes et des qualites- 
d'un systeme electrochimique (tel qu'une pile, une batterie, tin. 
systeme de dep6t elect rolytique, un systeme d 1 analyse'. 1 d'vii 
milieu) et la prevision de son fonctionnement futur, l ! uh des 
5 parametres dont la connaissance est essentielle est 1 r impedance- 
de ce systeme. T,"""*%. 

De fa<?on general e, un systeme electrochimique 
comprend, dans un milieu electrolytique, deux electrodes princi- 
pales - une electrode de travail et une cont re -electrode. Une 

10 electrode de reference est disposee au voisinage de 1* electrode 
de travail et est utilisee en relation avec celle-ci pour 
effectuer diverses operations de mesure ou de regulation. Dans 
la suite du present texte les tensions et courants mentionnes 
correspondent a des mesures effectuees entre une electrode de 

15 travail et 1* electrode de reference correspondante . 

L 1 impedance d'un systeme electrochimique correspond au 
rapport tension sur courant. II est eonnu que pour mesurer rapi- 
dement la valeur de cette impedance en fonction de la frequence, 
il est souhaitable de calculer le rapport des transformees de 

2 0 Fourier discretes (DFT) de la tension et du courant. Pour cela, 
on applique entre une electrode de reference et une electrode de 
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travail d'une cellule electrochimique un signal d» excitation de 

faible amplitude et on mesure a des intervalles de temps egaux 

AT, N valeurs successives e(n) et i(n) de la tension et du 

courant, avec 1 < n < N. L' expression generale des transformees 

de Fourier discretes E(f m ) et I(f m ) de la tension et du courant 

pour une sequence de N points est : 

N 

E(f m ) - AT ]T elnJexpC^TrjfmnAT) (D 



I(f m ) = AT 2 i(n)exp(-2njf m nAT) 



(2) 



ou j designe le nombre complexe dont le carre est -1, N le 
nombre de points de mesure, AT 1 1 interval le d ' echantillonnage . 
Le calcul de transformee de Fourier discrete peut etre effectue 
pour N/2 valeurs de frequence f m , avec 0 < m < (N/2)-l, N etant 
un nombre pair. Ces N/2 frequences discretes sont regulierement 
reparties entre 0 et 1/2 AT (0. . .1/mAT. . .1/2AT) . 

L' impedance complexe pour une frequence f m donnee est 
alors egale a : 

Z(f m ) = E(f m )/I(f ro ) - < 3 > 

Les methodes de mesure d' impedance utilisees dans la 
pratique, different principalement par le type de signal d 1 exci- 
tation du systeme : sinusoidal, multi-sinusoidal, bruit blanc, 
etc. Une excitation sinusoidale est de loin le procede le plus 
utilise car il s'avere etre le plus precis. Le procede consiste 
a imposer une regulation sinusoidale a la cellule electro- 
chimique avec un signal de petite amplitude et a enregistrer la 
reponse de la cellule en courant et en tension. Le rapport des 
transformees de Fourier de la tension et du courant a la 
frequence du sinus donnera la valeur de 1' impedance a cette 
frequence. Le spectre des frequences est balaye en modifiant la 
frequence d' excitation. 

Dans la theorie des systemes electrigues, 1' expression 
de l'inpedance (equation (3)) est correcte a condition que le 
systeme analyse soit lineaire et stationnaire et que des pertur- 
bations ne soient pas apportees par des phenomenes exterieurs au 
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systeme. Dans le cas des systernes electrochimiques , le respect 
de ces conditions impose des precautions particulieres . La 
fonction de transfert d'un systeme electrochimique est generale- 
ment non lineaire raais on peut la considerer conime lineaire sur 
5 une petite portion. C'est pour cela qu'on utilise un signal 
d» excitation de faible amplitude. En consequence, le rapport 
signal /bruit dirainue et on doit augmenter le temps de mesure, 
c'est-a-dire le norabre de points de mesure pour eliminer le 
bruit par une integration de la reponse. En revanche, 1' allonge- 

10 ment du temps de mesure pose des probl ernes quand a la stationna- 
rite du systeme Electrochimique. Dans beaucoup de cas reels, la 
condition de stationnarite ne peut pas etre remplie. Les causes 
sont multiples : relaxation du potentiel, relaxation du courant, 
relaxation des concentrations. Dans ces conditions, si temps 

15 de mesure est suf f isamment long pour que l'on ait un bori . rapport 
signal/bruit, le systeme n'est pas stationnaire pendant. Ma duree 
de mesure et l 1 impedance calculee n'a plus grande signification, 
notamment aux basses frequences (par exemple inferieures au 
hertz) ou les problemes de non- stationnarite sont particuliere- 

2 0 ment aigus. J 

Divers precedes ont ete proposes pour resoudre ce 
probleme de non stationnarite des systernes Electrochimiques, 
mais aucun n f a donne de solution satisfaisante. 

Ainsi, un objet de la presente invention est de 
25 prevoir un nouveau procede de calcul de 1' impedance d'un systeme 
electrochimique permettant de s'abstraire des erreurs liees a la 
non stationnarite du systeme. 

Pour atteindre cet objet, la presente invention prevoit 
un procede de determination de 1« impedance complexe Z(f m ) d'un 

3 0 systeme electrochimique non stationnaire, comprenant les etapes 

suivantes : 

mettre le systeme dans un etat de tension choisi et lui 
appliquer un signal sinusoidal de frequence f m , 
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mesurer immediatement apres des valeurs successives de 
la tension et du courant a des intervalles de temps reguliers 
AT, 

calculer les transformers de Fourier discretes de la 
tension (E(f)) et du courant (1(f)), la transformee de la 
tension etant calculee pour la seule frequence f m du signal 
sinusoidal et la transformee du courant etant calculee pour la 
frequence f m et pour deux frequences adjacentes f m _i et f ra+ i de 
part et d' autre de la frequence f m/ et 

calculer 1' impedance selon la formule suivante : 
Z(f m ) = E(f m )/I*(f m ) 
oil I* designe un courant corrige tel que : 

Re[I*(fm)] = Re[I(f m )] - {Re [I (f m+ i) ] +Re [I (f m -i> ] }/2 
lmtl*(f m )] - Im[I(f m )] - {lm[I(f m+1 )]+ImtI(f m -i)]}/2. 

Un mode de realisation particulier de la presente 
invention sera expose a titre non-limitatif en relation avec la 
figure jointe qui represente le spectre d» amplitude en fonction 
de la frequence d'un systeme electrochimique soumis a une 
excitation sinus oidale de faible amplitude et a un Echelon de 
tension . 

Tout d'abord, on notera que, conrme le potent iel 
applique est irrpose par 1" instrument d» analyse, la non station- 
narite du systeme ne peut se manifester que par des fluctuations 
de courant. 

La presente invention se base sur 1' analyse du compor- 
tement d 1 un systeme electrocnimique soumis a un echelon de 
tension. La DFT de la reponse en courant traduit povur 
l'essentiel des phenomenes de relaxation et done I'effet de la 
non stationnarite du systeme. Comme le represente la courbe en 
pointilles sur la figure unique, le spectre en amplitude I du 
module de la DFT du courant presente habitue llement une forte 
reponse a basse frequence a des echelons de tension. 

Par contre, la reponse a une excitation sinusoidale de 
frequence f m d f un systeme stationnaire se traduit par une raie 
unique a la frequence f m . En pratique, cette reponse correspond 
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d'une part a la reaction du systeme a 1* excitation a la 
frequence f m , d 1 autre part a la contribution des effets de 
relaxation. 

L'avantage de 1 'utilisation de la transformee de 
5 Fourier est que le spectre monof requence lie a I 1 excitation 
sinusoidale est superpose au spectre lie a 1' Echelon de tension. 

Selon la presente invention, une fois les considera- 
tions ci-dessus prises en compte, on propose de soustraire de la 
reponse en intensite a la frequence f m la contribution due aux 
10 relaxations, evaluee a partir de l f analyse de la DFT aux 
frequences adjacentes a f m . En effet, la transformee de Fourier 
pour des frequences f m -i et fm+i du systeme ne correspondra 
qu'aux phenomSnes de non stationnarit§ et on considerera que la 
valeur de l f amplitude liee aux non stationnarites pour une 
15 frequence f m est la moyenne des valeurs pour les deux frequences 
adjacentes. ,.: 
Ainsi, la presente invention propose : 

- de placer un systeme electrochimique dans un etat de 
tension choisi et de lui appliquer une excitation sinusoidale 

20 monof requence de faible amplitude, 

de mesurer, sans attendre la stabilisation du 
systeme, l r amplitude du courant et de la tension a des interval - 
les de temps reguliers immediatement apris 1 1 application de la 
tension tandis que I 1 excitation sinusoidale est appliquee, et 

25 - de calculer d'une part la transformee de Fourier 

discrete de la tension pour la frequence f m de 1' excitation et 
d 1 autre part la DFT du courant pour la valeur f m et pour deux 
valeurs de frequence f m -i et fm+i adjacentes a la frequence f m , 
et 

30 -de calculer la valeur de l 1 impedance complexe a la 

frequence f m a partir de la valeur de la TFD de la tension pour 
la frequence f ra et d'une valeur corrigee de la TFD du courant 
tenant compte des TFD calculees pour les frequences f m -i, f m et 
f m+l • 
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On calculera 1 • amplitude corrigee I* de la transformee 
de Fourier discrete du courant en valeur reelle et en valeur 
imaginaire par les deux equations suivantes : 

Re[I*(f m )] = ReEKf m )] - {Re [I (f m+1 ) ] +Re [I 1 }/2 

Hn[I*(f m )] = Im[I(f m )l - {iTtiEKfm+iJl+IraCKfm-i)]}^ 

La valeur de I 1 impedance corrigee des effets de non 
stationnarite est alors : 

Z(f m ) = E(f m )/I*(f m ) 
avec I*(f m ) = Re[I*(f m )] + jlm[l*(f m )] 

Des le releve des points effectue pour la frequence 
f m/ on pourra appliquer un signal sinusoidal a une nouvelle 
frequence et effectuer un nouveau releve, et ainsi de suite. 

On avantage du present procede est qu'il rend possible 
une analyse correcte d 1 impedance sur des systemes electro- 
chimiques non stationnaires notamment aux tres basses 
frequences. En meme temps, un gain de temps considerable est 
obtenu pour des systemes qui se stabilisent lentement car avec 
la correction selon l 1 invention, il n'est plus necessaire 
d'attendre la stabilisation apres une mise sous tension du 
systeme pour demarrer une analyse d'inp£dance. 
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FJSVSHDIC&TIONS 

1. Procede de determination de l 1 impedance complexe 
2(f m ) d'un systeme electrochimique non stationnaire, caracterise 
en ce qu'il comprend les etapes suivantes : 

mettre le systeme dans un etat de tension choisi et lui 
5 appliquer un signal sinusoidal de frequence f m , 

me surer immed i at eraent apr^s des valeurs successives de 
la tension et du courant a des intervalles de temps reguliers 
AT, 

calculer les transformees de Fourier discretes de la 
10 tension (E(f}) et du courant (1(f)), la transformee de la 
tension etant calculee pour la seule frequence f m du signal 
sinusoidal et la transf ormie du courant etant calculee ?'pour la 
frequence f m et pour deux frequences adjacentes f m -i et^f m+1 de 
part et d 1 autre de la fr§quence f ra , et £ 
15 calculer 1 T impedance selon la formule suivante 

Z(f m ) = E{f m )/I*(f m ) « 
ou I* designe un courant corrige tel que : T 
Re[I*(f m )] -Re[I(f m )] - {Re[I(f nH . 1 )]+Re[I(f m . 1 )]}/2 : - 
Im[I*(f m )] = lm[l(f ra )] - {llni[l{^ IH .i)]+lta[l(f m . 1 )]}/2? 
20 2* Procede selon la revendication 1, caracterise en ce 

qu'il est repete pour une succession de frequences d ' excitation . 
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